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1 . Thermal Instability in Loop Prominences 
(S.K. Antiochos and P.A. Sturrock! 

In our previous report we commented that in order for thermal 
instability to be important in loop prominences, it is necessary that 
radiation be the dominant cooling mechanism. In our model this is 
true because the magnetic field acts to decrease the conductive losses 
to the chromosphere. The reduction of conduction cooling is due to two 
effects: the curvature of the field increases the effective distance 

between the hot plasma and the chromosphere so that the heat flux is 
reduced; and the area of a loop increases reapidly with height so 
that for a large ratio of top area to base area, the flow of heat out 
of the ends of the loop is constricted. 

In order to investigate this effect quantitatively, we have 
considered a simple model In which the plasma remains stationary while 
cooling. Assuming that the temperature is high (~ 10 f K), gravity and 
radiation losses can be neglected. If s measures the distance along 
the loop measured from the top (Figure 1), the fluid equations reduce to 


n = n(s) 

(1) 

P = p(t) 

( 2 ) 

and the heat equation 


a -V/2 d p(t) _ 10 -6 d /A d \ . 

2 p dt A(s ) ds V d s P / * 

(31 

where 



(4) 
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For boundary conditions, wo assume that 


and 



( 5 ) 


T ( s ^ » t * = 'to 0 


( 6 ) 


d_T 

ds 


(O.t) = 0 


(7 * 


The Latter condition foLlows from the symmetry of the loop, since we 
expect the heat flux to vanish at the top. The variables s and t are 
separated in equation (3i, which can therefore be integrated directly 
to obtain 


P(t> “ P D C 1 + t/T) 


-2/5 


and 


T(X,t) = T(o ,o • (1 + t/T) 


•2/5 


1 - 


sin 2 V cos \ 


COS 


X COS Xjj 


Y 1 (x) 


(8) 

( 9 ) 


where the spatial variable has been changed from the distance s to the 
angle \ shown in Figure 1, The function as 1 for 0 5 X s T f/2, so 
that it can be neglected. The cooling time scale T can be calculated 
using the boundary conditions, 


T 


2i £o_^ sln2 x b 

10 10 -6 [t(o,o)} 7/2 cos U Xb 


( 10 ) 


To determine the effect of the area factor on the cooling, we 
have used the results above to compare the cooling time scale and 
outward heat flux of a loop of base area A to those of a loop of 
constant area but of exactly the same length. From equation (10) we 


obtain 


( 11 ) 


T /« of sin " *b 

T = \ s, / 4 

c \ b / cos y. 

where is the cooling time of a constant-area loop, and 

F_ I s b \ C03 *b 

F c \ H o h in \ 


(12) 


for the ratio of the outward heat fluxes. The right hand sides of 
equations (11 < and (12' can be expressed as functions of the base 
area only by using 


S b V \ 

H " 2 1 
o 


3v 2 + 1 


+ 


J&nfV'T V + l/jv 2 + 1) 

2^3 


(13) 


and 


where 


a - vh 3 

V^V 2 4- 1 


(14-. 


V = sin ^ 


(15) 


Figures 2 and 3 show the dependence of T/T c and F/? c on A^. We note 

that for small base areas the cooling time is greatly increased. 

Of course these results are not valid for very small areas since 

radiative cooling will begin to dominate conductive cooling. In order 

to accurately determine the heat losses, we ran several cases through 

our computer code (Antiochos and Sturrock, 1974 i. The results are shown 

in Figure 6. In each case the loop plasma was assumed to be initially 

7 

at rest with a uniform temperature and density profile of 10 K and 

11 -3 

10 cm respectively. The length of the loop from the top to the 
chromosphere was adjusted to be 5 * 10'* cm., while the values of A^ 
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that were used varied from I ./50 to L (corresponding to a loop of constant 
area) . 

From Figure b we notice that, as expected, loops of small base 
area have initially a lower cooling rate. However, after a large 
fraction of the loop energy has been lost the situation is reversed 
and loops of large base area cool more slowly. This is due to the 
following effect. At the later stages of cooling, the heat losses 
are essentially all radiative since conduction is negligible for low 
temperatures. Loops with small base area cool mainly by radiation 
and they tend to retain a uniform temperature and density profile as 
they cool. However, loops of large base area have significant conductive 
losses out of the base, and they tend to form large temperature differences 
between the feet and the top of the loop. Therefore, in order to 
maintain pressure balance the density of the hot plasma decreases, 
thereby lowering the cooling rate. 

This effect also occurs in thermal instability. The appearance 
of condensations does not necessarily mean that the loop cools faster. 

The c oling rate does increase at the onsei instability, but after 
a condensation has formed and absorbed a sigt.if j.cant fraction of the 
plasma, the hot gas remaining has a lower density and cooling rate. 

By incorporating this work with our previous results (Antiochos 
and Sturrock, 1975 1 on the thermal instability of loop plasma, we 
hope to explain the observed longevity of condensations in active 
prominences . 
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Figure Captions 

1. Geometry of loop prominences 

2. Ratio of conductive cooling rates t/t c for a dipole loop to that 
for a constant area loop of the same length as a function of A^, 
the ratio of base area to area at the top cv the dipole loop. 

3. Ratio of outward heat fluxes F/F as a function of A, . 

c b 

h. Fraction of energy left as a function of time for several loops 
with different base areas. 
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2 • F orce-Free Magnetic Fields and Their hole in Solar Activity 

. . W. Barnes ) 

During the present period, further calculations have been made 
using the spherical version of our force-free-field calculator 
discussed in the last report. Our model places a perfect conductor 
at the maximum radius of the mesh which prevents field lines from 
extending beyond that point. As long as the field strength there 
is small, this approximation does not introduce significant error. 

The cases reported last time showed a flattening of the field lineE 
against the boundary for large shear situations. We have increased 
the number of mesh points treated in the code so as to alleviate 
this problem. Also the nesh now extends from pole to pole, since 
it was found that field lines extend well into the southern 
hemisphere even though created near U 5 0 N latitude. The new mesh 
is 177 cells In £n r and 129 cells in 0. With the same length-to- 
width ratio of cell sides (0.7) as used before, the radius extends 
now to 20 surface radii instead of 5*2 radii. As before, the 
relaxation algorithm can be applied to the whole mesh or a coarser 
subset of it (covering the same physical region). The coarser meshes 
are used initially, with the finer meshes used for final convergence 
to the solution. 

The results cf a run which duplicates the surface conditions 
treated in the previous report are shown in Figures 1 - 4, The 
field lines are shown along with lines of constant shear (longitude), 
both radius and latitude are plotted on linear scales insteau of the 
radial logarithmic scale on which the calculation is performed. The 
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tt/ 2 , tt atuJ in radians. 


four figures show results lor shear angles of 0.*. 

The large degree to wb I tit the lines billow out with shear eatt be easily 
seen, it is also apparent that even with the outer boundary at £0 radii, 
the field lines for large shear eventually push against it. Figure 5 
shows the variation of stored energy with maximum shear angle y for 
these cases a td compares it with that of the open-field-line case. It 
can be seen that the mesh fineness causes some variation of stored 
energy due to inaccurate representation of the a function with coarse 
meshes. The solutions were ail self-consistent however, with calculated 
fields being very nearly force-free even with the coarsest mesh. 

We have also performed some calculations of fields in which a 
single field line is distorted outward at a single point. Figure 6 
shows the field pattern in the case with maximum shear = £tt. The 
local distortion of the or = 0,65 line can be seen. The result of 
this calculation unfortunately did not tell us whether the field energy 
increased or decreased with distortion, only that the change was small 
indeed (less than 1*$). Distorting the same field line in a case of 
small shear 0V max = TT /h), showed the same shallow energy dependence 
with a small increase in energy. Both of these cases had difficulty 
converging stably. We therefore feel that this approach will not be 
as useful as we had hoped. 

At present, we are trying a different approach to field line 
distortion. If the center of shear is placed at the equator and is 
made a line of symmetry, the equator becomes one of the model boundaries 
upon which we can arbitrarily set the value of the a- function. We 
feel that it will be easier to achieve convergence to a solution with 
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larger dLs Lor limits of the fields outward. If this is the case, we 
will be in a position to better understand the stability properties 
of cases with large shear angle. 


Figure C 
Figure l 


Figure 2 


Figure 3 


Figure 1+ 


Figure 5 


Figure 6 


ap Lions 

: Magnetic field lines (contours of Of) for a line dipole at 

b-5° N latitude on a sphere. The pattern is mapped onto 
a rectangle mesh with latitude plotted horizontally and 
radius plotted (logarithmically) vertically. The maximum 
shear in this case is small (0.1 radian). 

Magnetic field lines (contours of Of) and lines of constant 
shear (vi for the case - tt/2. 

Magnetic field lines (contours of of) and lines of constant 
shear (y ) for the case = tt. 

Magnetic field lines (contours of uf) and lines of constant 
shear (v) for the case = 2rt. 

Plot of relative stored magnetic energy versus maximum shear 
angle. Also shown is the energy associated with a dipole 
with completely open field lines. 

Magnetic field lines (df) and lines of constant shear (y) 

for the case v._ = 2 tt which has been artificially distorted. 

M 

The value of rt has been adjusted to represent the physical 
outward movement of a field line (that which shows a cusp). 
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Figure 5 
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3. Aircraft Accident Study 

(J.W. Knight and 1>,A. Sturrocki 

As we indicated in our previous progress report, we have nearly 
completed our study of the possibility that the aircraft accident rate 
may be affected by solar related phenomena. This study has so far 
produced provocative but not conclusive results. There are two extensions 
of our previous work which might make the results more clear-cut. 

Wilcox et al. (1973) found a dependence of the northern hemispheric 
vorticity area index on the interplanetary magnetic field structure. 

This effect was evident only for the winter months. If a connection 
between the aircraft accident rate and some solar related phenomena 
existed only for the winter months, our previous methods of analysis 
might not reveal it. Because the average accident rate is lower in 
the winter, any effect on the aircraft accident rate present only in 
the winter could be masked by the larger number of unaffected accidents 
in the summer months. To compute least-squares spectra for non-con- 
tinuous data sets, a major modification of our computer code which 
generates these spectra was necessary. We have recently completed 
this modification, and we intend to analyze the aircraft accident and 
Up data to ascertain if any seasonal effects are present. 

The earth's orbit is not circular, but slightly elliptical. This 
causes the apparent rotation period of the suri to vary through the 
year (Graf, 1974). The fact that the apparent rotation period of the 
sun varies through the year could make any solar-related signal present 
in the aircraft accident data more difficult to detect than a purely 
periodic signal. It is reasonably easy to derive, to first order in 
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e the eccentricity o£ the earth's orbit, a simple formula lor the 
amount cp the sun appears to have rotated as a function of time. If 
time is measured from perihelion, we find 


<*. = a> , . t ~ u) t - 2e sin to t = u>„ t - 24 sin uj t , 

’ cl n « n C^»n ' 


syn 


(U 


where w . , is the sideral angular velocity of the sun, is the average 

angular velocity of the earth in its orbit, and o# n Is the average 

apparent (synodic) angular velocity of the sun. It is relatively 

straightforward to modify the formula used to calculate the least 

squares spectrum to account for the eccentricity of the earth's orbit 

to first order in «. We simply vary the parameters a, b and cu to 

minimize 

N 

= N [x - <x> - a cos (cut - 2t sin oi t) - b sin (u)t - 2c sin u) t)]^ 
/ j n o o 

n=l 


V 3 


where x^ are the data, <x> is the average of the data, and N is the 
number of data analyzed. The significance of the least-square fits 
generated in this manner may be assessed in the same way as we have 
previously calculated significance estimates for least-squares spectra. 

We are in the process ot modifying the computer code used to calculate 
least-squares spectra to calculate least-squares "eccentricity corrected" 
spectra. 

We have used the computer codes written for the aircraft accident 
study to analyze the vorticity area index used by Wilcox et al. (1973) 
in their study of possible connections between the interplanetary 
magnetic field structure and terrestrial weather. Preliminary results 
of this analysis are discussed in Section 4 of this report. 
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h . Solar Effects on Terrestrial Weather 
(J.W. Knight and P.A. Sturrock'l 

It has long been known that geomagnetic disturbances tend to recur 
after approximately 27 days (Maunder, LyOha, lyjhb ; Chree and Stagg, 
1927; Bartels, 1932). More recently, this effect has been associated 
with the rotation of the interplanetary magnetic field structure (Ness 
and Wilcox, 1967; Schatten and Wilcox, 1067; Wilcox, I 968 ; Wilcox and 
Coburn, 1972). Evidence for some relation between the interplanetary 
magnetic field structure and terrestrial weather has also been reported 
(Wilcox et al., 1973; Wilcox et al., 1975). The interplanetary magnetic 
field is usually characterized by either two or four sectors, and this 
pattern rotates with a fairly well defined period of 27.1*. 1 days 
(Wilcox, 1968 ). The presence of a periodicity leads naturally to the 
use of spectral analysis as an independent test of the reported asso- 
ciation between the interplanetary magnetic field structure and terres- 
trial weather. We have obtained data for the vorticity area index 
used by Wilcox and his colleagues (Wilcox et al., 1973) in their 
analysis . 

Both the G and vorticity area indices exhibit annual variations 
P 

and long term trends. Both data sets have been normalized to remove 
these variations. We have calculated the least-squares spectrum of 
the vorticity area index and compared it with the least-squares spectrum 
of the 6 ^ index. The comparison was accomplished by tabulating the 
lesser of the two F statistics for the two spectra as a function of 
period, and searching for peaks in this "comparison spectrum". The 


five largest peaks in the comparison spectrum are displayed in Table 1, 



The largest peak occurred at a period of 27. ^9 days. We may use the F 
statistic to determine the probabilities that the calculated reductions 
in variance are due to chance. Following Abramowitz and Steguu (19703, 

we find the probability that the fit at 27*^9 days is due to chance for 

-b -c, 

the spectrum Is ~ 1.2 X 10 and for the vortictty is ~ 1.2 X 10 ' . 

The F statistic calculated from the data is distributed only approximately 

as M so using such low probability estimates is probably not 

justified. We therefore arbitrarily use the lowest probability tabulated 

- -k 

by Abramowitz and Stegun 1 1970 1 for the F distribution which is 10 , 

If the Gp and vorticity indices are independent, then the probability 

that both indices would exhibit fits of this significance at the same 

period is lQ Even if this e^t ‘mate is valid, it is not correct to 

use it since the period of 27.^9 days was not chosen a priori, but 

inferred from the data. It is more appropriate to assess the probability 

of finding one peak close to the "target" period of 27.1 days. Choosing 

the range in period to be 27*l ;±: .b days so that the observed peak is just 

included in the range is also improper. Rather, we choose the range of 

observed synodic rotation periods of the photosphere of the sun, 26.87 - 

28.9’j days (Newton and Nunn, I 95 I). The least-squares fits were 

calculated for 22 periods in this range. The probability of obtaining at 

-6 - c < 

least one of 22 probability estimates as small as 10 is ^ 2.2 X 10 
if all of the periods are independent. This is not the case, however, 
so the validity of the estimate is not clear-cut. 

The probability estimates calculated from the F statistic are 
unrealistically small and setting a lower limit of 10 ^ is not a very 
satisfactory procedure. We therefore have also calculated an empirical 
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probability estimate;. The F statistics for all l-T T periods analyzed 
were ranked, and the probability of each fit estimated to be the rank 
of the fit divided by the total number of periods for which fits wete 
calculated, To compare the spectra, the larger of these two probability 
estimates was tabulated as a function of period. The five most significant 
peaks in this comparison spectrum are displayed in Table 2. There were 
two periods with probabilities as low as or lower than the probability 
for the 27'h'9 day period. However, all three probabilities were close 
together and the next most significant peak had a probability estimate 
approximately twice as large. If we combine the two empirical probability 
estimates in the same manner as the estimates from the F test, we find 
that tne peak at 21. bQ days has the smallest probability estimate 
3 * 10 i of any of the periods analyzed bv more than a factor of 
two. Again we must compensate for the a posteriori choice of the period. 
The empirical probability estimate that the observed result was produced 
by chance is ^ 6.6 X 10 4 . 

The assumption that the G^ and vorticity indices are independent 
is probably not justified. However, it t-j difficult to estimate the 
probability that the calculated iits would appear in both spectra by 
chance without some Similar assumption. An alternate procedure is to 
assume that the relationship between the index and the interplanetary 
magnetic field structure is sufficiently well established that we may 
use the G^ spectrum to establish a rotation period. We may then accept 
the probability estimates for the 27.h- r * day period in the vorticity 
spectrum. These estimates are < IQ ^ utid ~ l.U X lcT 1 * for the F test 
and empirical probability estimates respectively. ft is possible that 

2h 


using the G^ spectrum to select a rotation period Is not appropriate* 

In this case, we may base our estimates only on the vorticlty spectrum, 
and the probability that the calculated tit would occur in one of 22 
periods is quite high (> for the empirical probability estimates. 


We have attempted to assess the significance of the observed peaks 
under various assumptions. The results vary from greater than 
to less than 75;i> We feel the former estimate is unrealistically high 
and the later unrealistically low. We feel the most reasonable procedure 
of the methods we have described is to accept the period derived from 
the Gp spectrum as solar related. We therefore feel that the results 
of our preliminary analysis are significant at approximately We 

do not feel that the preliminary analysis establishes a connection 
between the vorticity area index and solar activity; however, we do 
feel the analysis provides additional evidence for such a connection. 

We intend to check the results of our analysis by analyzing the 
vorticity data taking account of the eccentricity of the earth's orbit 
as discussed in Section 3 of this report. We also intend to explore 
the possibility of a the' theoretical methods of estimating the statistical 
significance of the results providing more realistic estimates. 



Table 1. 


Pe r i od 

F . 
min 

27.49 

9.07 

19.83 

8.37 

44.28 

7.58 

19.65 

6.87 

29.57 

5.77 


Table of the five largest peaks in a tabulation of the 

lesser of the two F statistics for both C and vorticity spectra. 

P 

Table 2 


Period 

Larger 
Probabi 1 ii 

19.88 

26/1921 

27.49 

27/1921 

44.52 

27/1921 

29.57 

53/1921 

37-58 

71/1^21 


Table of the five most significant peaks in a tabulation 
of the empirical probability estimates for both and vorticity 
spectra. 
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5. Stochastic Acceleration In Solar Flares 
(C. E. Newman ) 

The preliminary studies of the viability of stochastic acceleration 
by Mill) turbulence as a possible Phase 2 ion production mechanism in 
solar flares, as given in the last status report (Newman, 1975 )* have 
been modified to allow consideration of a particular model for the 
turbulent spectrum and supplemented by numerical calculations of the 
solutions. As mentioned in that report, a power- law frequency spectrum 
with index v for the turbulence results in an equation for the ion 
energy distribution which admits a power- law with index n - ~(v - l 1 
as a steady-state solution. The i *:tract lveness of stochastic acceleration 
as the ion production mechanism depends then upon how well it can provide 
answers to the following questions: 

1. For reasonable plasma parameters, characteristic of solar flares, 
does the ion distribution approach the power-law steady-state 
solution on a time-scale consistent with the observed 5-20 min. 
lifetimes of such bursts? 

2. Are the observed ion energies of 100 - 1000 toeV explained by this 
mechanism? 

3. Can it explain the enhanced abundances of heavier ions in these 
bursts? 

In order to obtain the answers to the above questions, we have 
assumed a specific model for the turbulence. In this model the frequency 
spectrum i' - a power- law with upper and lower cut-off frequencies cu. and 
respectively: 

. y/ 

Sio> “ u> for u>, - ui - at 
i u 
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S(cul - 0 for m < <u. or u> u> 

X u 


Since acceleration by M11D turbulence takes place via the quasl-llnear 
wave-particle interaction at cyclotron resonance (Newman, 1972), it 
follows that only ions with energies In the range E 6 E - ftu /ai.TE 
are accelerated by the turbulence. The iniection energy fin HeV) 
is given by 


E 

o 


2.5 X 10 


h 



t 


E 


f 


where Z and A are the charge and mass numbers of the ion, U - is 

H m p c 

the proton gyrofrequency , B is the magnetic field strength in gauss, 

-3 

and n g is the electron number density in cm 

The energy-space quasi-linear diffusion equation for the ion 
distribution N(E,t« becomes (Newman, 1973', under these assumptions 
and the additional condition that E^ > > , 


“ §£ N (E , t " = - (n + |l (E n N * + ^ (E n+1 N), E„ < E 


2 3E 


(13 


3E C 


where n = ~fv - 1), the diffusion coefficient D = E^ n t ^ 

2 n-1 o A 

and the characteristic time t^ (in seconds^ is given by 


t ^ o 
A 


x 10 


5 I “Kf £ ft) iRV 


(2> 


T 

The quantity — appearing in f2) is the ratio of the total energy 
*B b 2 

density of turbulence to the magnetic field energy density 7^ . 

The solution to (1) with boundary conditions 

p 

NfE, 0) = N (~1 , N(E , t) = N , N (°°, t) = 0 

O b o o 





was given in reference 1. This solution as a function of energy is 

plotted for various values of T = t/t in Figure 1, where n and cr 

A 

are taken as 2.5 aid 5*5 respectively. It is easily seen from this 
figure that the steady-state solution is nearly reached at time 
T = 0.1. 


It is now possible to provide the answer to the questions enumerated 
above. As parameters typical of a reconnection region, we take 


B n 30 gauss 


* ST " lo3 * “• 


9 „ -3 If _ 


10 cm 


ho . 


We find that the injection energy E q is thus given by 


E » 25 7 2 A -1 MeV , 
o 

and the highest energy attainable E^ by 

E, = 1000 Z 2 A" 1 MeV. 
t 


The characteristic time t. becomes 

A 

t A *. 5 X 10 3 Z A" 1 sec. 

so that the time required for a steady state to be achieved is 

_ i 

t ~ 0.1 t^ = 500 7 A sec. 

Thus the answers to questions (1) and (2) are definitely affirmative. 


To answer question (3)» we can show that for t & 0.1 t^ the 
total number density of particles accelerated above by this 
mechanism is given by 
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* 


iT - ^ t/t A + o.U , 

where 

N = / N (E,0) dE 

h ° 

o 

is the total number density of particles above E q at t = 0. Since 
varies as AZ ^ , we see that the heavier the ion, the greater the 
number (relative to the original density) accelerated above E q . 

Hence the last question can also be answered in the affirmative. 

For these reasons we conclude that MHD turbulence is an attractive 
mechanism for acceleration of energetic ions in Phase 2. However, it 
cannot explain acceleration of thermal ions (fa few KeV) ti. the 100 - 
1000 MeV range since the factor ( f/E^) in equation (2) becomes too 
large and gives values for t^ which are much too long. Thus some kind 
of injection (or preacceleration) mechanism is required to produce 
ions with energies ~ 25 MeV in order for stochastic acceleration to 
become effective. Feasible injection mechanisms will be evaluated 
and presented at a future date. 
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Figure Caption 

Figure l: The quantity N(E, t)/N o , where N is the solution to equation 

(l), is plotted versus E/E q for various values of t = t/t^. 
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6. Analysis of Hard X-Rays from Solar Flares 

B . J . Lipa and V. Petrosian 

In large solar flares, particles are accelerated in at least 
two distinct phases. In the flash phase, electrons are repeatedly 
accelerated over short time Intervals; In the second, they are 
accelerated to high energies, possibly by a shock wave which interacts 
with solar magnetic fields to create conditions under which Fermi 
acceleration occurs. Evidence for such a shock wave emerges from a 
study of pulsations which have been observed to modulate the radio 
continuum during solar flares (Achong, 197M- In the mechanism which 
has been proposed to explain these effects (Rosenberg, 1970), a Type II 
shock wave intercepts a magnetic flux tube that loops the corona. When 
the shock wave reaches the top of the arch, the tube is excited into 
radial oscillations and the resulting magnetic field oscillations 
produce pulsed modulation in synchrotron radiation. The same process 
can occur lower in the solar atmosphere, where modulation in the density 
produces pulsations in bremsstrahlung radiation. 

We have begun a study of hard x-ray emission, using data from 
the 0S0-5 satellite. To search for periodicities, maximum-entropy 
spectral analysis was used. Four flares have been analyzed and the 
results are summarized as follows. 

March 1 , I 969 

The flash phase consists of strictly periodic bursts of period 
39 seconds, which were observed to 111 KeV. Following the impulsive 
phase, thermal radiation dominates; no oscillations are observed. 
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March 21, 1969 

The flash phase consists of quasi-periodic impulsive bursts, 
separated by approximately 400 seconds. Superimposed on this were 
oscillations which were indicated by highly significant double peaks 
in the frequency spectrum with periods of 24 and 33 seconds. No 
oscillations were observed in a second acceleration phase. 

March 30, I 969 

The flash phase is extremely rapid and there is some evidence for 
short-period (8 seconds) oscillations. In the succeeding second 
acceleration process, no oscillations occur. 

November 1970 

Quasi-periodic bursts separated by approximately 14 q seconds 
dominate the first half of the event. Superimposed upon these are 

■f 

oscillations indicated by a highly significant peak in the frequency 
spectrum. The period varied between 29 and 48 seconds as the impulsive 
phase proceeded. 

It appears from this limited analysis that the first acceleration 
phase often consists of quasi-periodic bursts due perhaps to repetitive 
instabilities occurring in the magnetic field. Superimposed on these 
bursts are modulations which may be due to the modulation of the flux 
tube density described above. 

There is a difficulty in interpreting the periods found due to 
the manner in which the x-ray data was collected. The measurements 
were collected from the rotating wheel of the spacecraft by sampling 
the solar burst for ,19 seconds during every period of revolution of 
l.til seconds. It is therefore possible that the observed period is 
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a beat period with the revolution time. However, the periods obtained 
are of the same order of magnitude as the 16 -second period observed 
by Parks and Winckler (1971 K 

We intend to make a statistical study of x-ray data from a large 
number of flares. We intend to find the electron spectrum responsible 
for the x-ray bursts, the relation between the energy emitted and the 
periods involved, and the polarization of the resulting x-rays. We 
will also look for associated behavior in other solar phenomena. In 
particular we will look for pulsations at microwave and meter wave- 
lengths, and the relation between pulsating structures, Type II bursts 
and particle emission. 
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7. Solar X-Rays 
(V. Petrosian) 

At described in the previous semiannual report, the predictions of 
our model of the impulsive hard x-ray bursts are consistent with all 
observed properties of these bursts. This agreement is a consequence 
of the beaming of the x-ray-generating electrons toward the chromosphere. 
Because of the relativistic beaming of the bremsstrahlung radiation in 
the direction of electron motion, most of the radiation would be directed 
toward the chromosphere. This beaming becomes more important at higher 
energies. The reflection of the photons directed toward the sun could 
in certain cases contribute significantly to the observed radiation. 

We have devised a Monte Carlo procedure for calculation of the 
x-ray albedo of the solar atmosphere, for different x-ray energies and 
different incident angles. The accuracy of this procedure has been 
determined by comparing results from it with analytic solutions for 
a few cases where such solutions are possible. The results of this 
calculation would enable us to determine the center-to-limb variation 
with solar longitude of the energy spectral index, the frequency of 
occurrence, and the ratio of hard-to-soft x-ray intensities as 
described previously (Petrosian, 1975). These results are being 
prepared for publication. 

We intend to extend this program for calculation of polarization 
properties of the reflected radiation. There have been calculations 
of the polarization of x-rays from an electron beam (Haug, 1972 and 
Brown, 1972), but unfortunately only calculation of polarization of 
the photons emitted directly toward the earth has been published. 
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However, in order to determine the polarization state o£ the total (direct 
and reflected) x-rays we need to calculate first the polarization state 
of the x-rays emitted toward the sun. These calculations are under 
way. The results of this calculation will be used in the Monte Carlo 
program for determination of the polarization state of the reflected 
photons . 
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